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In  this  paper,  the  authors  investigate  a  robust  Integrated  Optimal  Design  (IOD)  devoted  to  a  passive  wind 
turbine  system  with  electrochemical  storage  bank:  this  stand  alone  system  is  dedicated  to  rural 
electrification.  The  aim  of  the  IOD  is  to  find  the  optimal  combination  and  sizing  among  a  set  of  system 
components  that  fulfils  system  requirements  with  the  lowest  system  Total  Cost  of  Ownership  (TCO).  The 
passive  wind  system  associated  with  the  storage  bank  interacts  with  wind  speed  and  load  cycles.  A  set  of 
small  power  passive  wind  turbines  spread  on  a  convenient  power  range  (2-16  kW)  are  obtained  through 
an  IOD  process  at  the  device  level  detailed  in  previous  papers.  The  system  cost  model  is  based  on  data 
sheets  for  the  wind  turbines  and  related  to  battery  cycles  for  the  storage  bank.  From  the  range  of  wind 
turbines,  a  “system  level”  optimization  problem  is  stated  and  solved  using  an  exhaustive  search.  The 
optimization  results  are  finally  exposed  and  discussed  through  a  sensitivity  analysis  in  order  to  extract 
the  most  robust  solution  versus  environmental  data  variations  among  a  set  of  good  solutions. 
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Nomenclature 

Neel  s 

number  of  cells  in  series 

Nch 

number  of  the  battery  bank  changes  over  20  years 

c 

scale  factor  (m/s) 

Ncyc 

equivalent  number  of  cycles 

Co 

cost  of  one  deep  battery  cycle  (€) 

Kyc 

equivalent  number  of  cycles 

C3 

battery  element  nominal  capacity  (Ah) 

Pext 

extracted  power  (W) 

Cbat 

battery  bank  cost  (€) 

Pload 

load  power  (W) 

(-Cell 

battery  cell  cost  (€) 

P\VT 

wind  turbine  nominal  power  (W) 

CF 

cycle  to  failure 

soc 

battery  cell  state  of  charge 

CwT 

wind  turbine  cost  (k€) 

Taut 

battery  bank  time  autonomy  (h) 

DOD 

depth  of  discharge 

TCO 

total  cost  of  ownership  (1<€) 

h 

battery  element  nominal  discharge  current  (A) 

V0 

battery  element  nominal  voltage  (V) 

I cel 

battery  cell  current  (A) 

V\X3 

generated  wind  speed  (m/s) 

Ich  max 

battery  cell  maximum  charge  current  (A) 

Vref 

reference  compact  wind  speed  (m/s) 

Idisch  max 

battery  cell  maximum  discharge  current  (A) 

Vwind 

wind  speed  (m/s) 

k 

shape  factor 

Wcyc 

weight  of  a  cycle 

N 

number  of  cycles  at  a  given  DOD 

T 

repeated  wind  cycle  (days) 

Neel _p 

number  of  cells  in  parallel 

Top 

operating  period  (days) 

1.  Introduction 

According  to  the  World  Energy  Agency  [1],  some  1.5  billion 
people  had  no  access  to  electricity  in  the  world  by  2009,  with 
more  than  80%  of  habitants  in  rural  zones.  Providing  consumers  in 
remote  areas  with  reliable  and  cheap  electricity  becomes  a  priority 
in  several  developed  and  undeveloped  countries  such  as  in  the 
case  of  isolated  cities  in  Tunisia.  The  steadily  increasing  demand  of 
fossil  fuels  along  with  concerns  about  global  warming,  presents 
natural  renewable  energy  sources  as  attractive  solutions.  Among 
these  sources  wind  energy  systems  (free  in  their  availability, 
renewable  and  non-polluting)  with  storage  are  among  the  most 
competitive  alternatives  for  electrifying  remote  consumers  and 
they  are  widely  used  in  both  autonomous  or  grid  connected 
applications.  These  systems  can  also  operate  in  parallel  with 
others  available  energy  sources  (fuel  cells,  diesel  generators)  and 
several  means  of  storage  (accumulators,  H2  storage,  etc.)  in  order 
to  enhance  the  system  reliability  [2-5], 

However,  the  drawbacks  of  such  sources  are  their  TCO  which  is 
still  expensive,  especially  for  small  wind  turbine  systems.  More¬ 
over  to  assure  the  service  continuity  and  to  protect  the  battery 
against  deep  discharges  (subsequently  extend  the  battery  bank 
life),  such  systems  require  an  additional  dynamic  source  of  energy 
or  an  optimal  wind  system  design  [8,9]).  Recently,  several 
researches  based  on  global  optimization  techniques  have  been 
focused  on  the  design  of  optimal  system  configurations  which 
meet  the  load  demand  for  given  climatic  data  [10-12], 

Bagul,  Borowy  and  Salameh  [6-8]  have  developed  several 
methodologies  for  optimally  sizing  a  wind/PV  system  associated 
with  a  battery  bank  for  a  given  load.  These  methods  are  based  on 
the  use  of  long  term  data  for  both  irradiance  and  wind  speed. 
However,  such  approaches  are  penalized  by  CPU  time  due  to  wide 
data  range.  Several  studies  have  used  the  average  hourly  wind 
speed  data  over  a  few  years  simulation  period,  but  this  vision 
strongly  filters  wind  turbine  powers.  Other  researchers  [13,14] 
have  developed  probabilistic  methods  to  determine  the  annual 
energy  of  a  wind  system.  In  particular  in  [15-19],  authors  have 
selected  the  optimal  combination  and  sizing  of  wind  generators, 
PV  modules  and  storage  batteries. 

This  paper  suggests  a  systemic  methodology  for  designing  the 
optimal  combination  and  sizing  of  passive  wind  turbines  asso¬ 
ciated  with  electrochemical  storage.  Generally,  deterministic  opti¬ 
mization  approaches  neglect  the  effects  of  environmental  inputs 
uncertainties  (including  variation  or  perturbation  of  wind  speed 
and/or  load  profiles)  which  can  lead  to  drastic  change  of  optimal 
solution  quality.  Several  studies  [20-26]  have  taken  into  account 


the  increased  need  for  sensitivity  analyses  to  perform  a  robust 
system  design  in  various  research  areas.Thus,  this  study  is  parti¬ 
cularly  focused  on  the  sensitivity  analysis  of  a  set  of  “good” 
solutions  to  obtain  an  optimum  system  quite  insensitive  to 
environmental  variations. 

The  remaining  of  the  paper  is  organized  as  follows.  The  passive 
wind  system  and  the  battery  characteristics  are  described  in 
Section  2.  In  the  third  section,  the  systemic  optimization  proce¬ 
dure  is  presented.  Section  4  is  dedicated  to  the  local  optimization 
procedure.  The  models  and  the  systemic  optimization  formulation 
are  exposed  in  Section  5.  Section  6  is  reserved  to  the  results  and 
the  sensitivity  analysis. 


2.  Description  of  the  system 

The  considered  system  is  a  “low  cost”  full  passive  wind  turbine 
(WT)  battery  charger  (Fig.  1)  without  active  control  and  with 
minimum  number  of  sensors  as  studied  in  [27,28]:  this  local 
optimization  loop  is  not  detailed  in  this  paper  but  only  referred 
to  previous  studies.  An  experimental  prototype  of  the  optimized 
system,  especially  the  PM  synchronous  generator,  has  been  built  in 
LAPLACE  Lab  and  had  confirmed  the  ability  of  the  passive  structure 
to  draw  a  power  close  to  the  optimal  range  with  small  losses.  This 
prototype  and  subsequent  study  is  detailed  in  [29],  The  wind 
turbine  parameters  have  been  obtained  by  applying  similitude 
relationships  with  reference  to  a  1.7  kW  wind  turbine  which  had 
been  previously  optimized  by  the  local  IOD  process  in  [30,31], 

Both  wind  speed  and  load  profiles  used  in  this  study  are 
considered  as  deterministic  data.  These  data  were  acquired  at  a 
typical  farm  in  Tunisia.  The  load  profile  is  set  on  24  h  and  day  by 
day  repeated  (Fig.  2).  The  wind  speed  profile  has  been  obtained 
from  a  previous  study  [32]  by  applying  a  “compact  synthesis 
process”  on  an  actual  wind  speed  profile  of  200  days  duration 
considered  as  reference  data,  with  the  aim  of  generating  a 


Wind  turbine  f  >( 


Fig.  1.  WT  system  with  battery  for  stand  alone  application  (rural  site 
electrification). 
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Fig.  3.  Wind  speed  compact  profile. 


Table  1 

Basic  characteristics  of  a  Yuasa  NP  38-121  lead  acid  battery  element. 


Nominal  capacity  C3  30.3  (Ah) 

Nominal  voltage  V0  12  (V) 

Nominal  discharge  current  J3  10.1  (A) 


compact  profile  on  a  reduced  duration  of  10  days  (Fig.  3)  to 
accelerate  the  optimization  process.  It  should  be  noted  that  this 
process  synthesizes  a  wind  cycle  which  fulfils  the  main  character¬ 
istics  of  the  reference  cycle  on  200  days,  especially  the  maximum 
wind  speed,  the  average  cubic  wind  speed  and  the  energy  content 
of  the  wind  cycle  [32]. 

in  this  study,  a  lead  acid  Yuasa  NP  38-121  [33]  is  considered  as 
battery  element.  The  basic  characteristics  are  summarized  in 
Table  1. 

NB:  in  the  following,  let  us  note  that  the  12  V  Yuasa  “battery 
element”  is  constituted  of  6  “cells”  of  2  V  in  series.  The  association 
of  several  “battery  elements”  (12  V)  in  series  and  in  parallel  will  be 
named  as  “battery  bank”. 

3.  Systemic  optimization  procedure 

In  order  to  perform  the  optimization  of  the  whole  system 
including  the  passive  WT  and  the  storage  battery  considering  wind 
speed  and  consumption  profiles,  we  have  adopted  an  approach 
based  on  two  optimization  levels: 

Level  1  Systemic  Optimization  (SO)  of  the  whole  system. 

Level  2  Local  Optimization  (LO)  of  the  WT  device. 

Given  a  wind  speed  profile,  the  SO  approach  consists  in 
selecting  a  WT  in  the  range  normally  provided  by  the  WT 
manufacturer  (WT  T]-Tn  and  the  corresponding  Permanent  Mag¬ 
net  Synchronous  Generator  (PMSG)  Gi-G„:  in  our  case,  these 
devices  are  obtained  from  the  LO  process  as  presented  in  [31]). 
The  couple  (T„  G,j  and  the  corresponding  storage  sizing  have  to 
satisfy  a  given  load  demand  at  lowest  TCO.  This  compromise  can 


WT  and  PMSG 


be  obtained  by  solving  the  optimization  problem  illustrated  in 
Fig.  4. 


4.  Local  optimization 

The  aim  of  this  second  optimization  level  is  to  build  a  range  of  n 
generators  corresponding  to  n  WTs  which  will  be  used  in  the  first 
optimization  level  (systemic  approach).  In  this  approach,  the 
battery  bank  voltage  is  as  constant  and  equal  to  48  V  (i.e.  a  series 
connection  of  4  Yuasa  12  V  elementary  batteries)  (Fig.  5). 

The  “mixed  reduced  model”  considered  in  this  optimization 
process  is  detailed  in  [27,31], 

In  [31  ],  an  IOD  method,  based  on  a  multiobjective  optimization, 
has  been  developed  for  sizing  the  elements  of  a  1.7  kW  passive 
wind  turbine  system.  The  WT  range  with  related  generator 
parameters  for  various  nominal  powers  has  been  obtained  by 
applying  similitude  relationships  with  reference  to  the  1.7  kW 
wind  turbine  system  [30],  Fig.  6  shows  the  extracted  powers  Pext  of 
the  passive  wind  turbines  (till  16  kW)  obtained  by  similitude  from 
the  optimized  reference  passive  wind  turbine  of  1.7  kW 
(Appendix).  It  can  be  seen  that  the  quality  of  wind  power 
extraction  of  these  passive  configurations  (solid  curves)  matches 
very  closely  the  behavior  of  active  wind  turbine  systems  operating 
at  optimal  wind  powers  by  using  an  MPPT  control  device  (i.e.  the 
cubic  curves  in  Fig.  6).  It  should  be  noted  that  the  selected  wind 
turbine  range  (1.7  up  to  16  kW)  is  constrained  by  the  average  load 
power  (Pfoad).  since  we  have  adopted  a  range  between  Vi  <Pload)  and 
5  < Pioad )•  Furthermore,  the  use  of  similitude  models  from  a 
reference  wind  turbine  device  of  1.7  kW  limits  the  range  of  power 
expansion  to  fulfill  validity  domain  of  similitude  relationships. 
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Fig.  6.  Extracted  power  of  the  range  of  passive  WT  systems  obtained  by  similitude  from  the  reference  1.7  kW  WT. 


5.  Systemic  optimization 

The  aim  of  the  SO  stage  is  the  minimization  of  the  TCO  on  a  life 
cycle  of  20  years  of  the  passive  WT  coupled  with  a  storage  bank 
ensuring  the  electrification  of  the  isolated  farm  under  a  specific 
compact  wind  speed  cycle  [32].  To  achieve  this  optimization 
process,  a  cost  model  of  each  system  component  has  to  be 
determined. 

5.1.  WT  cost  model 

Generally,  the  WT  subsystem  is  composed  of  a  turbine,  a 
nacelle,  a  tower,  a  rectifier,  a  mechanical  transmission  system. 
There  is  no  single  component  that  dominates  the  WT  cost.  Typical 
owning  costs  given  by  “eaglewestwind"  [34]  for  a  range  of  turbines 
from  2  kW  up  to  20  kW  are  shown  in  Fig.  7:  system  costs  include 
costs  due  to  all  components  of  the  device.  A  linear  interpolation  of 
these  costs  has  led  to  the  following  WT  cost  model: 

CWt  =  1-7Pwt+3  (1) 

The  lifetime  of  WTs  being  assumed  to  be  at  least  20  year,  the 
owning  cost  of  WTs  will  simply  be  due  to  the  investment  cost, 
which  assumes  that  no  maintenance  costs  are  needed  due  to  the 
robustness  of  the  full  passive  wind  turbine  structure.  Such 
assumption  can  be  justified  for  low  power  wind  turbines,  espe¬ 
cially  as  the  passive  structure  is  particularly  robust. 


Fig.  7.  Assumption  of  a  WT  cost  model  versus  rated  power  Pwj. 


5.2.  Battery  bank  cost  model  based  on  cycling  effects 

The  battery  cost  depends  on  the  charge  and  discharge  cycling 
in  terms  of  number  and  cycle  depth  that  the  battery  can  provide 
(battery  lifetime  strongly  depend  on  both  depth  and  rate  of 
discharge).  The  equivalent  number  of  cycles  (Ncyc)  is  obtained  by 
a  weighted  sum  of  the  cycles  in  terms  of  their  depth  of  discharge 
(DOD).  The  weights  are  derived  from  the  manufacturer  curve 
giving  the  number  of  complete  cycles  that  the  storage  element  is 
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capable  of  delivering  in  terms  of  its  DOD  (characteristic  called 
“cycle  to  failure"). 

The  lifetime  model,  for  Yuasa  NP  38-121  lead  acid  battery,  uses 
an  exponential  curve  fitting  based  on  the  available  cycle  to  failure 
(CF)  curve  (Fig.  8)  versus  DOD  [35,36,37]: 

CF  =  1 77.77  +  7807.39.e-6  75DOD  (2) 


Considering  the  number  of  cycles  provided  by  the  battery  with 
a  DOD  of  100%  as  a  reference,  the  weight  of  a  cycle  depth  DOD 
(■ wcyc(DOD ))  is  expressed  by  the  following  equation: 


Wcyc(DOD)  = 


CF(100%) 
Cf(DOD ) 


(3) 


Thus,  the  equivalent  number  of  cycles  is  obtained  by: 

Ncyc  =  wCyC(DOD)  x  N(DOD)  (4) 

DOD 


where  N(DOD)  is  the  number  of  cycles  at  a  given  DOD  obtained  by 
the  Rainflow  method. 

For  an  operating  period  rop  of  20  years  and  a  repeated  wind 
cycle  t  of  10  days,  the  approximate  cost  of  the  battery  bank,  taking 
account  of  cycling  constraints  over  20  years,  is  expressed  by  the 
following  equation: 

CBAT(k€)  =  Ncel  sNcei  pC0  x  1 0 -3NlyCT-f  (5) 

where  Ncei_s  is  the  number  of  battery  cells  of  2  V  associated  in 
series  (here  considered  as  constant  and  equal  to  24  to  constitute  a 
battery  of  48  V),  Ncei_p  is  the  number  of  battery  cells  associated  in 
parallel  and  C0  is  the  cost  of  one  deep  battery  cycle  (DOD=  100%) 
estimated  at  0.1  €  [33].  This  latter  cost  comes  from  the  Yuasa  12  V 
battery  element  cost  which  is  108  €.  Then,  the  2  V  battery  cell  cost 
is  CCeii=  108/6  €  for  180  allocated  deep  cycles.  Finally,  one  cell  deep 
cycle  cost  is  CCe///180=0.1  €. 

It  should  also  be  noted  that  the  initial  purchase  cost  of  the 
battery  CBato  is  defined  as  following: 

Cbato(M€)  =  Ncei_sNcei_pCceu  x  10  3  (6) 
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Fig.  9.  Systemic  process  of  optimal  design. 


compact  profile  are  then  extrapolated  over  20  years  in  order  to 
estimate  the  TCO  due  to  cycling  effects. 


5.3.1.  TCO 

A  single  objective  function  is  used  consisting  in  the  TCO  which 
includes  WT  and  battery  bank  owning  costs  over  a  duration  period 
of  20  years. 

TCO  =  Cwt  +  Cbat  (7) 


5.3.  Optimization  problem  formulation 

The  problem  is  to  develop  a  systemic  approach  that  allows 
optimizing  system  configurations  (passive  WT  with  storage  bank) 
that  satisfy  load  power  demand  with  minimum  TCO.  The  TCO 
calculation  includes  the  WT  owning  cost  and  battery  bank  owning 
cost.  The  systemic  optimization  process  is  detailed  in  Fig.  9.  By 
means  of  the  compacted  wind  cycle  obtained  from  the  synthesis 
process  proposed  in  [32]  (as  presented  in  Fig.  3),  only  10  days  of 
system  operation  have  to  be  simulated  in  this  optimization 
process.  Results  obtained  from  this  typical  but  representative 


Depth  of  discharge  DOD 

Fig.  8.  Cycle  to  failure  curve  for  Yuasa  accumulator  cell. 


5.3.2.  Design  variables 

The  optimization  problem  only  uses  two  discrete  design  variables: 

•  N cei_p :  The  number  of  battery  cells  associated  in  parallel  (NceJ_s_ 
being  the  number  of  battery  cells  (2  V)  in  series  is  fixed  at  24 
for  a  48  V  battery  bank). 

•  index  i:  This  index  identifies  the  couples  of  turbines  and  associated 
generators  (7)  and  G,j  used  in  the  simulation  model  for  computing 
the  system  TCO  and  the  optimization  constraints.  Note  also  that 
this  index  corresponds  with  power  level  (in  kW)  of  the  WT. 

5.3.3.  Optimization  constraints 

•  g, :  is  a  constraint  related  to  the  maximum  discharge  current 

( ldisch_  max)- 

The  battery  cell  maximum  discharge  current  (max(/cei)),  must 
be  less  than  the  current  Idisch_max- 
g,  =  ma  x(Icei)-Idisch 

max  —  0  (8) 

•  g2:  is  a  constraint  related  to  maximum  charge  current  ( ICh_max )■ 
The  absolute  value  of  the  maximum  batteiy  cell  charge  current 
( |  min(/ce;)  | )  must  be  less  than  the  current  ICh_max'- 

g2  =  min (Icei)-lch  max  <  0  (9) 

Note  that  charge/discharge  current  are  respectively  considered 
as  negative/positive. 
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•  g3:  is  a  constraint  related  to  the  battery  cell  state  of  charge 

(SOC). 

The  minimum  value  of  the  battery  cell  state  of  charge  min(SOC 
(t))  must  be  greater  than  20%: 

g3  =  20%—  min(SOC(t))  <  0  (10) 

In  this  study,  Idisch_max  and  lCh_max  have  been  chosen  equal  to  the 
nominal  discharge  current  /3  =  10.1  A  (Table  1). 

As  defined,  the  optimization  problem  is  sufficiently  simplified 
(one  objective,  two  design  variables  and  three  constraints)  to 
justify  the  use  of  an  exhaustive  search  (80.14=  1120  combinations) 
instead  of  adopting  more  sophisticated  optimization  algorithms. 
Other  approaches  may  be  used  in  more  complex  case  study,  as  for 
hybrid  systems  with  several  sources.  Note  also  that  this  simplified 
and  efficient  optimization  problem  is  mainly  due  to  the  structure 
of  the  design  process  with  2  loops,  one  for  the  WT  device  and 
another  one  for  systemic  approach. 

6.  Results  and  sensitivity  analysis 

6.1.  Optimization  results 

Fig.  10  shows  a  set  of  solutions  resulting  from  the  systemic 
optimization  process,  corresponding  to  a  TCO  lower  than  115.5  k€. 
It  should  be  observed  that  the  TCO  range  is  displayed  with  a 
reduced  scale  (1/115  k€)  so  that  TCOs  of  these  solutions  can  be 
considered  as  quasi-equivalent.  Each  point  corresponds  to  one 
solution,  (i.e.  one  couple  of  variables:  index  i  representing  the  WT 
rated  power  and  Nce(_p).  The  cheapest  solution  (circled  in  Fig.  10) 
corresponds  to  a  configuration  with  63  branches  in  parallel  (i.e. 
252  “12  V  Yuasa  battery  elements”  for  the  battery  bank)  with  a  WT 
of  13  kW  of  nominal  power  (i=13).  Table  2  gives  the  different 
characteristics  of  some  quasi-equivalent  optimal  solutions  and 
Fig.  11  shows  the  state  of  charge  SOC  and  the  battery  cell  current 
(/ce/)  evolutions  of  the  solution  Sol  13opt.  One  can  see  that  a 
reduced  DOD  range  is  obtained  from  that  optimal  sizing  taking 
account  of  cycling  effect  on  system  costs. 


Note  that  NCh  is  the  number  of  the  battery  bank  changes 
necessary  for  20  years.  NCh  can  be  calculated  from  (5)  and  (6)  by 
deriving  the  following  equation: 

M  CMT  over  20  years 

Nch  =  - r -  (11) 

i-BATO 

Another  criterion  necessary  for  the  optimal  solutions  analysis  is 
the  battery  bank  autonomy  time  Taut  of  each  solution.  This  time  is 
defined  as  the  corresponding  duration  of  a  complete  battery 
discharge  up  to  20%  of  SOC,  considering  the  load  profile  without 
any  wind  power  production.  Table  3  displays  the  autonomy  time  of 
each  solution. 

Fig.  12  confirms  the  example  of  Fig.  11  by  illustrating  a  major 
trend:  optimization  of  TCO  strongly  limits  variations  of  battery 
SOC  (reduced  values  of  DOD)  in  order  to  enhance  its  lifetime. 
Results  would  have  been  completely  different  if  the  cycling  was 
not  taken  into  account  through  the  system  TCO.  In  such  a  case,  the 
whole  range  of  DOD  would  be  exploited  with  deep  discharge 
cycles. 

The  analysis  of  the  results  summarized  in  Table  2  illustrates 
that  some  optimized  solutions  with  different  characteristics  have 
almost  the  same  TCO.  The  solutions  indexed  15(a),  15(b) 
(Pwt=15  kW)  and  16(a),  16(b)  (PWT=  16  kW)  correspond  with  the 
cheapest  battery  owning  cost  CBat  among  others  configurations. 
The  choice  of  the  solution  indexed  13  is  optimal  in  terms  of  TCO, 
but  nearly  equivalent  and  surely  not  the  most  relevant  solution 
regarding  other  criteria  as  the  system  autonomy.  Fig.  12  shows  the 
curves  of  the  probability  density  functions  {pdf}  of  the  number  of 
cycles  versus  the  depth  of  discharge  (DOD)  and  their  locations  on 
the  corresponding  “cycle  to  failure”  curve  for  each  solution.  The 
circle  diameter  on  the  “cycles  to  failure”  curve  is  proportional  to 
the  cycle  occurrence,  i.e.  the  equivalent  pdf.  These  figures  show 
that  the  maximum  of  DOD  for  all  solutions  does  not  exceed  23%  of 
the  total  storage  capacity  of  the  battery  bank.  These  results  explain 
why  there  is  a  trend  to  increase  the  number  of  battery  cells  in  the 
solution  Sol  15(a)  compared  with  the  solution  Sol  16(a):  this  is 
obviously  to  further  expand  the  total  storage  capacity  of  the 
battery  bank  and  to  reduce  subsequently  deep  discharge  cycles 
therefore  enhancing  the  lifetime. 


Table  2 

Quasi-equivalent  optimal  solutions  characteristics. 


Sol 

16(a) 

Sol 

16(b) 

Sol 

12 

Sol 

13(a) 

Sol 

130pt 

Sol 

15(a) 

Sol 

15(b) 

Neel _p 

49 

59 

61 

62 

63 

69 

79 

Cwt  (1<€) 

30 

30 

23 

25 

25 

28.5 

28.5 

Nrvr  (on  20  vears) 

714 

594 

622 

602 

589 

520 

452 

Cbat  (k€) 

84.1 

84.2 

91 

89.7 

89.1 

86.2 

85.8 

TCO  (k€) 

114.3 

114.4 

114.5 

114.8 

114.2 

114.7 

114.3 

CftATO  (k€) 

21 

25.4 

26.3 

26.7 

27.2 

29 

34.1 

NCh  (on  20  years) 

3 

3 

3 

3 

3 

2 

2 

100 

90 


u 

S  80 


v\rv  \ 

riwy 

\A 

TO 

IN 

V 
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Fig.  11.  SOC  and  Icei  evolutions  for  the  optimal  solution. 


Table  3 

Time  autonomy  of  quasi-equivalent  optimal  solutions. 


Sol  16(a)  Sol  16(b)  Sol  12  Sol  13(a)  Sol  13opt  Sol  15(a)  Sol  15(b) 

Taut  (h)  5  6  6.2  6.4  6.6  7.1  8.1 


Fig.  10.  Similar  solutions  with  equivalent  TCO. 
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Table  4 

Wind  speed  profiles  characteristics. 


c 

k 

Vmin 

V max 

<V> 

<v3> 

Vref 

9.5 

2.3 

0.1 

25 

8.3 

1195 

Vl 

9 

2 

0.04 

25.3 

8 

867 

v2 

10 

2 

0.04 

31.6 

8.9 

1189 

v3 

12 

2 

0.1 

34.5 

10.8 

2072 

Furthermore,  at  equivalent  system  cost,  it  is  certainly  prefer¬ 
able  to  dispose  of  a  greater  storage  capacity  in  order  to  ensure  load 
demands  during  longer  periods  without  wind  (i.e.  autonomy 
Table  3).  More  generally,  it  is  interesting  to  analyze  these  solutions 
by  sensitivity  analysis  versus  changes  of  environmental  data  (wind 
speed  and  load)  then  to  rebuild  the  optimization  problem  to 
choose  a  robust  and  optimal  solution. 

6.2.  Sensitivity  analysis  of  optimal  solutions  versus  environmental 
data  variations 

This  section  is  devoted  to  the  sensitivity  analysis  of  the 
equivalent  optimal  solutions  presented  in  Fig.  10  versus  environ¬ 
mental  input  changes  (i.e.  wind  speed  and  load  profile  variations) 
with  the  aim  of  identifying  the  most  robust  solutions. 

6.2.3.  Sensitivity  versus  wind  speed  profile  changes 

Equivalent  optimal  solutions  represented  in  Fig.  10  have  been 
obtained  for  a  particular  wind  speed  compact  profile  (Vrt,f)  (see 


V(m/s) 

Fig.  13.  Probability  distribution  functions  of  the  wind  speed  profiles. 

Fig.  3).  Vref  distribution  corresponds  to  the  characteristics  of  a 
Weibull  law  with  a  scale  factor  c=9.5m/s  and  a  shape  factor 
k=  2.3.  Here,  we  aim  at  analyzing  solution  robustness  versus 
statistic  characteristics  mentioned  in  Table  4. 

For  this  purpose,  three  particular  time  profiles  of  wind  speed  with 
200  days  duration  are  generated  from  different  Weibull  distributions 
(V],  V2 ,  V3)  as  displayed  in  Fig.  13  and  characterized  in  Table  4. 

Table  5  gives  the  number  of  cycles  over  20  years  NCyo  the 
battery  bank  cost  CBat  and  the  TCO  of  several  equivalent  optimal 
solutions  for  each  wind  speed  profile.  These  results  show  that  the 
solution  Sol  15(b)  is  the  most  robust  in  terms  of  TCO  change:  (Sol 
15(b)  remains  lower  than  the  one  obtained  with  other  solutions 
whatever  the  wind  conditions.  This  robustness  is  partly  due  to  the 
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Table  5 

Wind  speed  profile  variations  results. 


Sol  16(a) 

Sol  16(b) 

Sol  12 

Sol  13(a) 

Sol  13opt 

Sol  15(a) 

Sol  15(b) 

Kef 

Ncyc 

714 

594 

622 

602.98 

589 

520 

452 

CBat  (k€) 

84.1 

84.2 

91.1 

89.7 

89 

86.2 

85.8 

TCO  (k€) 

114.32 

114.40 

114.38 

114.86 

114.25 

114.74 

114.37 

V , 

Ncyc 

1015 

829 

995 

924 

900 

728 

623 

CBat  (k€) 

119 

117 

145 

137 

136 

120 

118 

TCO  (k€) 

149 

147 

169 

162 

161 

149 

146 

v2 

Ncyc 

846 

659 

806 

752 

736 

582 

480 

Cbat  (k€) 

99 

93 

118 

111 

111 

96 

91 

TCO  (k€) 

129 

123 

141 

137 

136 

124 

119 

v3 

Ncyc 

653 

455 

557 

499 

492 

402 

340 

Cbat  (k€) 

76 

64 

81 

74 

74 

66 

64 

TCO  (k€) 

107 

94 

104 

99 

99 

95 

92 

Table  6 

Wind  speed  profile  V  characteristics. 


c 

k 

^ min 

^ max 

<V> 

<v3> 

V 

9.5 

2.3 

0.18 

25.67 

8.41 

914 

Days 


Fig.  14.  Load  profile  changes. 

battery  capacity  of  the  Sol  15  (b)  (NCei_p= 79)  which  is  the  highest 
among  all  other  solutions. 

6.2.2.  Sensitivity  versus  load  profile  changes 

The  sensitivity  analysis  of  the  optimal  solutions  versus  load 
profile  is  performed  by  considering  an  actual  wind  profile  of  one 
year  duration  (V)  whose  characteristics  are  given  in  Table  6. 
Changes  made  with  respect  to  the  original  load  profile  consist  in 
inserting  a  variation  of  25%  of  load  power  during  the  even  days 
and  -25%  during  odd  days  (Fig.  14). 

Table  7  shows  that  the  Sol  15(b)  remains  the  most  robust 
solution  in  terms  of  TCO  versus  variations  of  load  profile.  As 
mentioned  in  the  case  of  wind  speed  profile  variations,  this  is 
explained  by  the  highest  storage  capacity. 

Beyond  this  sensitivity  analysis,  we  can  guess  that  the  solution 
with  the  highest  battery  capacity  (Sol  15  (b))  is  also  the  most 
autonomous  in  case  of  wind  drops:  it  is  subsequently  the  most 
reliable  and  relevant  for  this  application. 

7.  Conclusion 

This  paper  illustrates  a  systemic  optimization  approach 
devoted  to  the  optimal  design  of  a  full  passive  WT  with  storage 


Table  7 

Sensitivity  versus  load  profile  variations. 


Sol  16 

(a) 

Sol  16 

(b) 

Sol 

12 

Sol  13 

(a) 

Sol 

130pt 

Sol  15 

(a) 

Sol 

15(b) 

Original 

Ncyc 

824 

655 

842 

761 

746 

578 

495 

profile 

Cbat 

(1<€) 

97 

92 

123 

113 

112 

95 

94 

TCO 

(k€) 

127 

122 

146 

138,4 

138 

124 

122 

Modified 

Ncyc 

791 

650 

811 

746 

728 

573 

479 

profile 

Cbat 

(k€) 

93 

92 

118 

111 

110 

94 

90 

TCO 

(k€) 

123 

122 

142 

136 

135 

123 

119 

dedicated  to  stand  alone  applications  (here  for  rural  electrification 
purpose).  System  environment,  especially  wind  and  load  condi¬ 
tions  is  integrated  in  the  systemic  design  by  optimization.  The 
integrated  optimal  design  problem  was  divided  into  two  optimiza¬ 
tion  processes:  a  local  optimization  which  aims  at  designing  a  set 
of  optimized  WT  in  order  to  constitute  a  manufacturer  range  to  be 
used  in  the  second  optimization  process  dedicated  to  the  whole 
system  design.  The  systemic  optimization  objective  is  to  minimize 
the  total  owning  system  cost,  integrating  simultaneously  WT  and 
battery  costs  by  taking  account  of  the  number  of  device  change 
due  to  cycling  effects.  It  is  interesting  to  note  that  such  a  systemic 
vision  of  the  TCO  over  the  whole  life  cycle  (here  20  years)  leads  to 
oversize  batteries  in  order  to  reduce  DODs  (lower  than  20%)  in 
order  to  reduce  the  number  of  device  changes  over  the  20  years. 
Optimal  results  were  analyzed  and  discussed  then  completed  by  a 
sensitivity  analysis  which  has  compared  several  selected  solutions 
previously  considered  as  quasi-equivalent.  The  sensitivity  of  these 
quasi-optimal  solutions  has  been  analyzed  versus  changes  of 
environmental  data  (wind  speed  and  load  profiles  variations) 
showing  that  solutions  with  highest  storage  capacity  were  the 
most  relevant,  being  the  most  robust  and  also  the  most  autono¬ 
mous  in  case  of  lack  of  wind.  The  availability  of  the  developed 
methodology  was  successfully  demonstrated  with  the  considered 
environmental  data  and  it  is  appropriated  for  other  acquired  sets 
of  experiment  data  to  find  consistently  the  optimal  generation 
system  design. 
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Table  Al 


Turbine:  rated  power 

Radius  (m) 

Friction  coefficient  (N  m  s/rd) 

Reference:  1.7  kW 

1.25 

0.025 

T,:  1 

0.95 

0.0087 

T2:  2 

1.35 

0.034 

T3:  3 

1.66 

0.07 

T4:  4 

1.91 

0.13 

T5:  5 

2.14 

0.21 

T6:  6 

2.34 

0.31 

T7:  7 

2.53 

0.42 

T8:  8 

2.71 

0.55 

T9:  9 

2.87 

0.70 

T,0:  10 

3.03 

0.86 

T„:  11 

3.17 

1.04 

T,2:  12 

3.32 

1.24 

T,3:  13 

3.45 

1.46 

T,4:  14 

3.58 

1.69 

T,5:  15 

3.71 

1.94 

T16:  16 

3.83 

2.21 

Table  A2 

PMSG:  rated  power  (kW) 

Flux  (Wb) 

Resistance  (m£l)  Inductance  (mH) 

Reference:  1.7 

0.20 

12626 

1.34 

C,: 

1 

0.15 

278.4 

1.57 

C2: 

2 

0.21 

94.3 

1.36 

C3: 

3 

0.28 

60.9 

1.53 

G4: 

4 

0.33 

40.4 

1.5 

C5: 

5 

0.35 

26.5 

1.33 

Gq’. 

6 

0.42 

24.2 

1.55 

C7: 

7 

0.41 

15.6 

1.22 

C8: 

8 

0.47 

14.6 

1.37 

G9: 

9 

0.53 

13.7 

1.51. 

G10 

10 

0.47 

8.3 

1.05. 

C„ 

11 

0.51 

7.9 

1.14. 

G,  2 

12 

0.56 

7.6 

1.22. 

Gi3 

13 

0.61 

7.3 

1.3 

G14 

14 

0.65 

7 

1.38 

C,5 

15 

0.7 

6.8 

1.46 

C,6 

16 

0.75 

6.6 

1.54 

Appendix 

See  Tables  Al  and  A2. 
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